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ABSTRACT 

The development of a system to detect the radio-frequency (RF) pulse as- 
sociated with extensive air showers of cosmic rays is described. This work 
was performed at the CASA/MIA array in Utah, with the intention of de- 
signing equipment that can be used in conjunction with the Auger Giant 
Array. A small subset of data (less than 40 out of a total of 600 hours of 
running time), taken under low- noise conditions, permitted upper limits to 
be placed on the rate and strength of pulses accompanying showers with 
energies around 10 17 eV. 

1 Motivation 

As a result of work in the 1960's and 1970's [TJ E3 E| , some of which continued beyond 
then (see, e.g., [HI E] ) , it appears that air showers of energy 10 17 eV are accompanied 
by radio-frequency (RF) pulses [9], whose properties suggest that they are due mainly 
to the separation of positive and negative charges of the shower in the Earth's magnetic 
field |10[ ITT] . The most convincing data were accumulated in the 30-100 MHz frequency 
range. However, opinions differ regarding the strength of the pulses, and atmospheric and 
ionospheric effects have led to irreproducibility of results. In particular, there may also 
be pulses associated with cosmic-ray- induced atmospheric discharges [121 US] • There are 
reports of detection at MHz or sub-MHz frequencies (HI El El, which could be associated 
with such a mechanism. Signals above 100 MHz have also been reported [5|. 

A study was undertaken of the feasibility of equipping the Auger array Jl] with the 
ability to detect such pulses. The higher energy of the showers to which the array would 
be sensitive could change the parameters of detection. Before a design for large-scale 
RF pulse detection could be produced, it was necessary to retrace some of the steps of 
the past 30 years by searching for the pulses accompanying 10 17 eV showers, and by 
studying some of the factors which led to their irreproducibility in the past. RF pulses 
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may be able to provide auxiliary information about primary composition and shower 
height PJ. 

In this article we describe the prototype activity at the CASA/MIA site and draw 
some conclusions from it regarding plans for the Auger project. We have not been able 
to demonstrate the presence of RF pulses at CASA/MIA, and could only set upper limits 
for their intensity. The upper limits on the rates of events in which the North-South (or 
East- West) projection of the signal pulse was greater than some value were established 
at R up (S uNS > €°) = 0.555/(£°) 2 h" 1 and R up (S uEW > S°) = 0.889/(£°) 2 h" 1 , 
respectively, with £ v being the electric field strength per unit of frequency, measured in 
^V/m/MHz. More concrete plans for RF detection at Auger must await a prototype 
at the Auger site which utilizes some of the lessons learned from the present work. A 
preliminary description of this work was presented in Ref . ^5] • 

In Sections 2 and 3 we discuss expectations for RF signals and previous claims of 
their observation. Section 4 is devoted to details of the setup at CASA/MIA, including 
some of the reasons for choosing the specific configuration utilized. Our preliminary 
results are given in Section 5, while Section 6 deals with issues specific to a giant array 
such as Auger. We summarize in Section 7. Appendix A gives details of the sensitivity 
calculation, Appendix B establishes some properties of simulated pulses, while Appendix 
C summarizes cost estimates for an installation at the Auger site. 

2 Expectations 

We briefly summarize some expectations [3] for the characteristic of the RF signal asso- 
ciated with cosmic ray air showers. 

2.1 Mechanisms of pulse generation 

In the 1950's, R. R. Wilson [TJ] proposed that cosmic rays could induce the atmosphere 
to act as a giant spark chamber, triggering discharges of the ambient field gradient. This 
gradient, normally around 100 V/m, can attain values as high as 10 kV/m during intense 
thunderstorm activity [3]. Thus, the mechanism would lead to pulses of greatly varying 
intensity, whose correlation with air showers would be difficult to establish unless field 
gradients could be monitored over the whole path of the discharge. 

Another mechanism of pulse generation is associated with the asymmetry in electron 
and positron yields in showers as a result of Compton and knock-on processes. By the 
end of the shower, electrons outnumber positrons by about 10-25%, leading to a transient 
of vertically moving negative charge jU]. This is thought to be the main mechanism for 
generation of radio-frequency signals from showers in solid material such as polar ice 
[To] or sand ^7j, but is probably not the dominant mechanism in the atmosphere. 

Still another source of electromagnetic radiation in a cosmic ray shower involves 
the separation of positive and negative electric charges in the Earth's magnetic field. 
This is thought to be the dominant mechanism accounting for atmospheric pulses with 
frequencies in the 30 - 100 MHz range, and will be taken as the model for the signal for 
which the search was undertaken. 
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2.2 Characteristic pulse 



The time profile of a pulse due to charge-separation in the Earth's magnetic field can 
be modelled |S 3 by assuming that the bulk of the shower giving rise to the pulse is 
concentrated between an altitude of 10 and 5 km (for a shower of energy 10 17 eV) and 
calculating the pulse duration by comparing the total path lengths between the antenna 
and the beginning and the end of the shower. The rise time of a pulse from a shower 
with zero zenith angle observed 200 meters from its axis is expected to be about 5 ns, 
followed by a longer decay time and a still longer recovery time with opposite amplitude 
(about 100 ns) such that the total DC component is zero. 

The radiation from any stage of the shower which is traveling directly toward the 
antenna is expected to arrive to the antenna about the same time as the shower itself. 
The difference is accounted for by the refraction index of air. Such an essentially 5- 
function pulse has the highest-frequency components in its spectrum. Showers for which 
the impact parameters of the cores are farther from the antenna will have reduced high- 
frequency components since the total pulse duration will be longer, approaching several 
microseconds for vertically incident showers viewed from the side. 

The pulse is expected to grow linearly in amplitude with shower energy as a result 
of the increased number of particles emitting RF energy. This linear growth assumes 
coherence of the emitting particles, which is probably a good assumption for RF wave- 
lengths of several meters. The greater penetration of the atmosphere by more energetic 
showers also leads to an increased RF signal since the radiating particles are closer to 
the receiver. This should make the pulse amplitude increase more rapidly than linearly 
with primary energy. However, this effect is largely offset by the fact that a greater frac- 
tion of such deeply penetrating showers will have reduced high-frequency components 
in their pulses, as a result of the greater apparent time taken by the pulse to build up 
to its maximum amplitude at the receiver. The combination of the above three effects 
is expected to lead fortuitously to an overall linear dependence of pulse amplitude on 
primary energy 

At extremely high energies, shower particles will even be lost by collision with the 
Earth. This may give rise to a different type of RF signal but will not be effective in 
the context of the charge-separation mechanism considered here. The pulses associated 
with charge separation in the Earth's magnetic field should correspond to radio signals 
with approximately horizontal polarization. (For showers not arriving vertically from 
directly overhead there will also be a small vertical polarization component.) 

2.3 RF backgrounds 

Discharges of atmospheric electricity constitute an important source of background 
pulses. These will be detected at random intervals at a rate which depends strongly 
on local weather conditions as well as on ionospheric reflections. Man-made RF sources 
include television and radio stations, police and other communications services, broad- 
band sources (such as ignition noise), and sources within the experiment itself. (We shall 
discuss such sources for the CAS A/MIA array presently.) The propagation of distant 
noise sources to the receiver is a strong function of frequency. During years of sunspot 
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minima (e.g., 1995-6), ionospheric propagation on frequencies above 25 MHz is rare 
except for "sporadic-E" propagation, which can permit signals to arrive from distances 
of up to 2000 km via a single reflection from the ionosphere. As solar activity increases 
(e.g., subsequently to 1996), consistent daytime propagation over even greater distances 
can occur on frequencies up to and beyond 30 MHz. 

Galactic noise can be the dominant signal in exceptionally radio-quiet environments 
for frequencies in the low VHF (30-100 MHz) range PJ. For higher frequencies in such 
environments, thermal receiver noise becomes the dominant effect. We shall see that the 
CASA/MIA site is far from quiet enough that these effects become limiting. 

3 Some previous observations 

An early proposal involved detection of the ionization produced by air showers via radar 
[T%1 119j. The first claim for detection of the charge-separation mechanism utilized rela- 
tively narrow-band techniques at 44 and 70 MHz [T|. A Soviet group reported signals at 
30 MHz [20j, while a University of Michigan group at the BASJE Cosmic Ray Station on 
Mt. Chacaltaya, Bolivia j2] studied pulses in the 40-90 MHz range. The collaboration of 
H. R. Allan at Haverah Park in England [3] studied the dependence of signals on primary 
energy E p , perpendicular distance R of closest approach of the shower core, zenith angle 
6, and angle a between the shower axis and the magnetic field vector. Their results 
indicate that the electric field strength per unit of frequency, S v , could be expressed as 



where Rq is an increasing function of 9, equal (for example) to (110 ± 10) m for v = 
55 MHz and 6 < 35°. 

The Haverah Park observations are consistent with the model mentioned in Sec- 
tion 2.2 in which the pulse's onset is generated by the start of the shower at an elevation 
of about 10 km above sea level, while its end is associated with the greater total path 
length (shower + signal propagation distance) associated with the shower's absorption 
about 5 km above sea level. (The elevation at the CASA/MIA site is about 1450 m 
above sea level; the average atmospheric depth is 870 g/cm 2 |2T|.) 

The Haverah Park observations were subsequently updated to yield field strengths 
approximately 12 times weaker than Eq. (JTJ [4j, while observations in the U.S.S.R. gave 
field strengths approximately 2.2 times weaker than ([!}. Thus, some question persists 
about the magnitude of the effect, serving as an impetus to further measurements if the 
RF detection technique is to be employed as part of a new giant array. 

More recent pulse detections include claims for pulses with components below 500 kHz 
seen by observers at the AGASA array in Akeno, Japan [5] and a group working at 
Yakutsk in Siberia [7j, and claims for pulses at VHF frequencies seen by groups at the 
Gran Sasso in Italy jH] and at Gauhati University in India [H] . There seems to be no una- 
nimity regarding the time duration, generation mechanism, or intensity of these pulses. 
Related methods have been used to study lightning-induced pulses 
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Figure 1: Geometry of the CAS A/MIA array. Small squares denote CASA stations; 
cross-hatched rectangles denote muon patches. Large rectangle near center is the central 
trailer; rectangle to right (east) of the array is RF trailer; symbol "A" denotes placement 
of antenna. 

4 CAS A/MI A Prototype setup 

4.1 Description of the CASA/MIA detector 

The Chicago Air Shower Array (CASA) [2T| was originally constituted as a rectangular 
grid of 33 x 33 stations on the surface of the desert at Dugway Proving Ground, Dugway, 
Utah. The inter-station spacing is 15 m. A station has four 61 cm x 61 cm x 1.27 cm 
sheets of plastic scintillator each viewed by its own photomultiplier tube (PMT). When 
a signal appears on 3 of 4 PMTs in a station, a "trigger request pulse" of 5 mA with 5 /xs 
duration is sent to a central trailer, where a decision is made on whether to interrogate 
all stations for a possible event. Details of this trigger are described in Ref. |21j . 

For future reference, we shall denote the coordinates of each box by (n x ,n y ), where 
— 16 < (n x ,n y ) < 16, n x = (x/15 m), n y = (y/15 m), and (x,y) denotes the position 
of the center of the box to the (East, North) of the center of the array. When this 
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experiment was begun the CASA array had already been reconfigured to remove boxes 
with —16 < n x < —13, i.e., the 4 westernmost "ribs" of the array. For runs performed 
in 1998, boxes with n x = 16 had also been removed from the array. 

The University of Michigan collaborators designed and built a muon detection array 
(MIA) to operate in conjunction with CASA. It consists of sixteen "patches," each 
having 64 muon counters, buried 3 m below ground at various locations in the CASA 
array. (See Fig. 1.) Each counter has lateral dimensions 1.9 m x 1.3 m. Four of the 
patches, each about 45 m from the center of the array, lie on the corners of a skewed 
rectangle; four, each about 110 m from the center of the array, lie on a quadrangle with 
slightly different skewed orientation, and eight lie on the sides and corners of a rectangle 
with sides x ~ ±180 m and y ~ ±185 m. 

In April of 1991 the CASA/MIA array was partially disabled by a lightning strike 
which hit one of the few trees on the site. The array was repaired, and an extensive 
lightning-protection grid installed. The grid consisted of wires strung on poles about 
15 feet above the array, traveling in the x, y, and x ± y directions. This grid turned out 
to have significant effect on our choice of parameters for the RF studies. 

During the operation of the present experiment, 144 surface Cerenkov detectors |2*3] 
were distributed throughout the array. Other additions to the array, which shall not 
concern us, included a stereographic atmospheric Cerenkov detector (DICE) [21] and 
an optical facility for communicating with a high-resolution atmospheric fluorescence 
detector (HiRes) located on a hilltop several miles away. 

4.2 Initial RF surveys at CASA/MIA site 

In order to determine whether RF pulse detection was feasible at the CASA site, a 
spectrum analyzer was used to make a broad survey of the RF noise at the CASA 
site in various frequency ranges and at various locations. It was determined that in the 
central trailer, the broad-band noise associated with various computers, switching power 
supplies, and other electronics was so intense that no RF searches could be undertaken. 
The same was true to a great extent at any position within the perimeter of the lightning- 
protection grid. Moreover, it was deemed unsafe to erect an antenna above that grid 
within the perimeter of the array, since any projecting object would defeat the purpose 
of the grid. 

Surveys just outside the array indicated a much quieter RF environment. An antenna 
was placed about 24 m east of box (16,0), corresponding to x = 263.8 m, y — m, and 
its signal fed into a trailer located about 10 m closer to the array. All further studies 
were performed using this configuration. (See Fig. 1.) Nonetheless, there still remained 
a number of identifiable noise sources, which we now describe. 

4.2.1 Television and FM broadcast stations 

The CASA/MIA site is about 100 km southwest of Salt Lake City, at first sight affording 
a reasonably quiet RF environment. However, many television and FM stations in Salt 
Lake City broadcast from a high mountain about 35 km southwest of Salt Lake City, 
or 65 km northeast of Dugway. These are responsible for a major component of the RF 
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Table 1: Television stations broadcasting from site 65 km northeast of Dugway. 



( Q 1 1 Gl rrn 


wildlliltJl 


JjcLIIU. 


V IUxaJ 




1 UWcI 






(MHz) 


(MHz) 


(MHz) 


(kW) 


KUTV 


2 


54-60 
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KTVX 
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KSL-TV 
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166.0 


KBYU-TV 
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198-204 


199.25 


203.75 


162.0 


KSTU 


13 


210-216 


211.25 


215.75 


112.0 



signal in the range which is of greatest interest to us. As an example, we summarize 
the VHF television stations broadcasting from the above site in Table 1. The video 
and audio frequencies shown are carrier frequencies. Video signals are modulated with 
vestigial-sideband modulation, occupying the range from 1.25 MHz below the carrier 
frequency to about 3.5 MHz above it. A color subcarrier lies 3.58 MHz above the video 
carrier. Audio signals are frequency modulated with deviation not exceeding 250 kHz 
so as to remain within the total allotted bandwidth of 6 MHz for each channel. The FM 
broadcast band, extending from 88 to 108 MHz, is packed with strong signals, with the 
strongest typically spaced by the 0.8 MHz interval characteristic of inter-station spacing 
in a large urban area. 

4.2.2 CASA noise 

The CASA boards contain crystals oscillating at various frequencies, including 16, 20, 
and 50 MHz. The behavior of a single CASA board was investigated at the University 
of Chicago. The various clock signals were detected at short distances (< 1 m) from 
the board, but a much more intense set of harmonics of 78 kHz emanated from the 
switching power supplies. These harmonics persisted well above 100 MHz. At 144- 
148 MHz (monitored using an amateur radio transceiver), they overlapped, leading to 
intense broad-band noise. 

The above signals were considerably less problematic at the RF trailer. During CASA 
operation the boards' clock frequencies and some of their harmonics (including 32, 40, 
and 48 MHz) were detectable. However, noise from the switching power supplies seemed 
to be at an acceptably low level. 

The CASA boards emit powerful RF pulses when digitizing and transmitting data. 
These pulses constituted a major background to our RF search, and will be discussed in 
Section 5. The noise arrived through the antenna system and not through the trigger 
cable or antenna feed cable, as was determined by acquiring data with a dummy load 
in place of the antenna. 
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4.2.3 Intermittent narrow-band interference 

In addition to persistent RF carriers from TV and FM broadcast stations, intermittent 
signals would appear from time to time. The strongest of these was traced to local 
narrow-band FM communications. This signal was so strong that digital filtering meth- 
ods (to be described below) were powerless to eliminate it. Consequently, any event 
containing such a signal was discarded for further analysis. 

4.2.4 Low-frequency interference sources 

Although the majority of survey work dealt with frequencies above 25 MHz, some effort 
was made to reproduce claims of low- frequency ("LF") pulses [S], which for our purposes 
will be taken to involve frequencies below 500 kHz. (The AM broadcast band contains 
numerous signals above 530 kHz, preventing the study of higher frequencies.) Initial 
surveys were performed using a Sony SW-7600G all-band portable receiver and an ICOM 
IC-706 amateur transceiver. However, considerably greater sensitivity was achieved 
using a Palomar VLF converter which converts the band 10-500 kHz to 3510-4000 kHZ, 
which was then detected using the IC-706. 

The major source of interference at the site was a nondirectional aircraft beacon 
(NDB) operating at the Dugway airport on 284 kHz. Other NDBs and other LF carriers 
above about 110 kHz were detectable but considerably weaker. A custom-made filter 
was procured j2Z] to suppress the carrier at 284 kHz and signals from the AM broadcast 
band above 500 kHz. This filter was employed during some of the low-frequency studies 
to be described below. 

4.3 Measurement considerations and initial setup 

As mentioned above, the location of the receiving antenna (about 30 m east of the edge 
of the CASA array, at x = 263.8 m, y — m), was dictated by a compromise between 
proximity to the array and reduction of noise. This noise was carried, to a large extent, 
by the lightning-protection grid which overlays the array. 

It was decided at an early stage to concentrate on the search for horizontally polar- 
ized pulses as described in Section 2. Consequently, a broad-band antenna with linear 
polarization was adopted. Initial surveys were taken with one of the original antennas 
from the Mt. Chacaltaya experiments j2], which had been preserved from the 1960s. 
This antenna was a large model manufactured for VHF television reception, with some 
elements which had been added by the experimenters to improve low-frequency response. 

The Mt. Chacaltaya antenna was mounted on a portable searchlight tower attached 
to a small trailer. The tower could be extended to a height of about 35 feet. The 
antenna was slightly damaged in a collapse of the tower as a result of improper latching 
procedures. As insurance against further such incidents, a portable military surplus 
log-periodic antenna was acquired. This antenna (a Dome and Margolin model to be 
described below) was found to have superior response in the frequency range of interest 
and very robust construction (even surviving a subsequent tower collapse), and was 
adopted for subsequent studies. 
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Table 2: Modes of filtering, (a) Suppression at 284 kHz and above 500 kHz in some 
runs. 



Mode 



3 dB bandpass (MHz) 



Narrow-band 
Broad-band 
Low-frequency 



25-35 
25-250 
0.05-2.5 (a) 



The antenna was mounted on the fully-extended searchlight tower with its center 
at a height of 35 feet above ground, with the favored direction of reception arriving 
from the zenith, and with arbitrary azimuthal orientation. Data were taken with two 
orientations: "East- West" polarization and "North-South" polarization, both referred to 
magnetic North (14° east of true North [21] at Dugway). In addition, a projecting arm 
of the mounting bracket was used to suspend a 10-meter-long vertical antenna which 
was used for the low-frequency surveys. 

The bandwidth to be covered by the RF search was not initially specified, but to be 
determined by experience with survey experiments. Consequently, two main modes were 
used, a narrow-band mode and a broad-band mode. These are compared in Table 2, 
where we also list a low-frequency mode used in the LF survey. Their implementation 
is described in Sec. 14.4.31 

Some previous investigations (e.g., j2] and |3|) were able to detect RF pulses using 
a "stand-alone" trigger based on the reception of transients alone. This possibility 
was investigated using a broad-band receiver with filters admitting several different 
frequency ranges, and demanding coincidences of signals received in a minimum number 
of channels. It was found that the vast majority of such transients at the Dugway site 
were not associated with CASA/MIA events; they were probably due to atmospheric 
discharges. Such "stand-alone" transients, in fact, were found to increase during periods 
of enhanced atmospheric electrical activity. As a result, our main results concern RF 
data taken with a trigger based on large CASA/MIA events. We comment further on 
the possibility of a "stand-alone" trigger for future experiments in Section 6.4. 

The trigger was formed at the central CASA/MIA trailer, in a manner to be described 
in detail below. It was communicated to the RF trailer over RG-59 cable. The electrical 
length of the cable was found to correspond to a pulse delay of 2.15 fis. No evidence 
for pickup of this trigger pulse from the antenna was found. Other methods considered, 
and rejected in favor of the simpler electrical communication, included optical fiber and 
infrared sensors. 

4.4 Design features 
4.4.1 Antenna system 

A portable log-periodic antenna manufactured by Dome and Margolin, Model # DM 
ARM 160-5, with a nominal response of 30-76 MHz, was acquired from FairRadio Co. in 
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Table 3: Properties of log-periodic antenna used for RF studies 



Nominal frequency range (MHz) 30-76 

Usable frequency range (MHz) 28-170 
Number of elements 9 

Dimensions (m) 3x3 

Feedline RG-58U 60 feet 



Table 4: Properties of preamplifiers used for RF studies 



Manufacturer 


Model 


DC power 


Gain 


Frequency 






(V) 


(dB) 


range (MHz) 


Mini-Circuits 


ZFL-500LN 


13.6 


26 


DC-500 


ANZAC 


AM-107 


18 


10 


1-500 



Lima, Ohio, for about $60. (A spare was used for noise studies at the University of 
Washington.) Overload protection was provided by two 1N4148 diodes of opposite po- 
larity connected across the antenna terminals, leading to a maximum output voltage of 
about ±0.6 V. A gas discharge tube manufactured by Alpha/Delta provided lightning 
protection. Some properties of the antenna are summarized in Table 3. 

4.4.2 RF front-end 

The RF amplification stage consisted primarily of one or two ZFL-500LN low-noise 
broad-band preamplifiers manufactured by Mini-Circuits, and for certain runs low-noise 
preamplifiers manufactured by ANZAC. Specifications of these preamplifiers are sum- 
marized in Table 4. 

4.4.3 Filtering 

Table 5 contains a summary of all filters used in the experiment with the exception of the 
284 kHz filter [27| described previously. These filters are manufactured by Mini- Circuits; 
they were obtained with tubular cases fitted with BNC connectors. 

A typical "narrow-band" configuration described in Table 2 involved feeding the 
signal from the antenna through the feedline, a BHP-25 filter and a BLP-30 filter with 
combined 3 dB points of 25 and 35 MHz, a ZFL-500LN preamplifier with 26 dB of 
gain, a BBP-30 filter with 3 dB points 25 and 35 MHz, another Mini-Circuits ZFL- 
500LN preamplifier with 26 dB of gain, and a BHP-250 filter to suppress any high- 
frequency noise. (Some data runs involved permutations of these components. The 
above configuration was found to minimize feed-through of preamplifier noise. Some 
runs involved a dual ANZAC preamplifier instead of a ZFL-500LN.) 

A "broad-band" configuration involved the same feedline and BHP-25 filter, a sin- 
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Table 5: Filters used in RF data acquisition 



Model 


Type 


3 dB point (s) 






(MHz) 


BLP-1.9 


Low-pass 


2.5 


BHP-25 


High-pass 


25 


BLP-30 


Low-pass 


35 


BBP-30 


Bandpass 


25, 35 


BLP-250 


Low-pass 


250 



gle ZFL-500LN preamplifier, and a BLP-250 filter. A "low-frequency" configuration 
involved the feedline, a BLP-1.9 filter and a ZFL-500LN preamplifier, with a 284 kHz 
notch filter inserted before or after the BLP-1.9 in some runs. The notch filter also 
contained a roll-off above 500 kHz. 

4.4.4 "Large-event" trigger and design 

A trigger based on the coincidence of seven of the eight outer muon "patches" (see 
Fig. 1) was set to select "large" showers in the following manner. Each muon patch 
was set to produce a trigger pulse of length 5 /xs and amplitude —120 mV when n 
of its 64 counters registered a minimum-ionizing pulse within 5.2 fis of one another. 
For engineering runs (until 12/28/96), n was set equal to 4, while for later runs it was 
increased to 5 to favor larger showers and reduce noise. The pulses were then combined 
in two groups of 4, feeding through two 2X attenuators into two fan-in/fan-outs (to avoid 
saturation of inputs) and the resulting pulses further combined to produce a summed 
pulse. This signal was fed to a LeCroy 821 Discriminator, whose output was amplified 
to an amplitude of about —6 V and then sent over RG-59 cable to the RF trailer (see 
Fig. 1). The trigger pulse at the RF trailer had an amplitude of about —2.4 V and a 
duration of 1 /is. 

The above trigger was estimated to correspond to a minimum shower energy of 
somewhat below 10 16 eV, based on the integral rate [2H] at 10 18 eV of 0.17/km 2 /day/sr. 
At this level good correlation could be established between trigger pulses and events 
recorded by the CASA data acquisition system. 

Only shower radiation that is stronger than 3 /iV/m/MHz can exceed the average 
noise level at CASA site by three standard deviations or more [Sec. 4.6]. According to 
the original Haverah Park results [Eq. (JTJ)], a typical shower that would lead to such 
radiation is a vertical shower of energy 10 17 eV or higher at a distance of 210 m. If the 
rate for showers with energy greater than E behaves as 1/E 2 , showers above 10 17 eV 
would be expected to occur with a rate of 17/km 2 /day/sr. Since RF detection relies on 
muon triggering, the antenna can only detect radiation from those showers whose cores 
pass inside the rectangle of the array. Just about 0.06 km 2 of the array area lies inside 
the 210 m radius from the antenna. With the solid angle of observation limited by a 
zenith angle of 50°, one expects about 2.25 detectable shower pulses per day. 
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4.4.5 Data acquisition 

A Tektronix TDS-540B digitizing oscilloscope registered filtered and preamplified RF 
data on a rolling basis. These data were captured and stored on hard disk using a Na- 
tional Instruments GPIB interface upon receipt of a large-event trigger. Data were taken 
using various computers at different times, allowing for analysis both at the University 
of Washington and at Chicago. The Washington system used a Macintosh Quadra 950 
running Labview, with a latency time of about 8 seconds between events, while the 
Chicago system used either a Dell XPS200s desktop or a Dell Latitude LM laptop run- 
ning a C program adapted from those provided by National Instruments, with a latency 
time of about 2 seconds. Each trigger caused 50 /is of RF data, centered around the 
trigger and acquired at 1 GSa/s, to be saved. 

4.4.6 Rates and off-line processing 

The total trigger rate ranged between about 20 and 50 events per hour, depending on 
the value of n = 4 or 5 of muon counters chosen to generate a patch trigger pulse and 
on intermittent sources of noise sometimes present in the trigger system. Concurrently, 
the CASA on-line data acquisition system was instructed via a special program called 
MUTRIG to write files of events in which at least 7 out of the 8 outermost muon patches 
produced a patch pulse. These files, one for each CASA run, typically overlapped 
with the records taken at the RF trailer to a good but not perfect extent [3U] as a 
result of occasional noise on the trigger line. Moreover, an undiagnosed timing problem 
occasionally caused the loss of a muon trigger pulse for certain large events recorded by 
MUTRIG. 

4.4.7 Off-line overload rejection 

Events were typically recorded at a gain such that the maxima and minima corresponded 
to about 2/3 of the dynamic range of the oscilloscope's 8-bit data acquisition system 
(ranging from —128 to 127 digitization units). Local intermittent monochromatic RF 
signals occasionally saturated this dynamic range. Such events were rejected off-line by 
discarding any cases in which the maxima and minima exceeded 100 digitization units. 
In the configuration used for the final bounds on pulse amplitudes, corresponding to an 
oscilloscope setting of 5 mV per division of 25 digitization units, we thus rejected all 
signals corresponding to preamplifier peak outputs greater than ±20 mV. In some cases, 
with oscilloscope settings of 20 mV per division, we rejected signals with preamplifier 
peak outputs greater than ±80 mV. In all cases these voltages were well below the 
manufacturer's specified limit of ±400a/2 mV (3 dBm), and within a satisfactorily linear 
range of preamplifier response. 

4.4.8 Calibration 

The average gain G ant of the antenna in its forward direction rises from about 3 dBi 
(decibels with respect to an isotropic radiator) at 30 MHz to a peak of about 5 dBi at 
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50 MHz, slowly decreasing to 4 dBi at 76 MHz |3T]. We shall take an average gain of 
4 dBi (G ant = 2.5) over the frequency range of interest. 

More precise calibration would involve modelling of the gain pattern using a program 
such as EZNEC [32]; an d integrating over directions of expected signal arrival. This 
modelling also would have been useful in order to emulate the frequency-dependent phase 
distortion induced by the antenna, but it was not found possible to obtain a sufficiently 
close fit to the antenna's measured standing-wave-ratio characteristics to take such a 
model seriously. Certainly this point should be addressed in any future studies. One 
might also utilize sources of known strength such as amateur radio satellites broadcasting 
on 29.4 MHz, FM and television stations, and galactic and solar noise. The existing 
data contain signals from FM and television stations broadcasting near Dugway, some 
of whose field strengths are well enough known that they may be usable for calibration. 
Alternatively, for future work it would be helpful to calibrate antennas on an antenna 
range at some distance from an impulse generator, broadcasting through a broad-band 
antenna with already-determined characteristics. 

4.5 Signal processing 
4.5.1 Fourier methods 

In order to remove strong Fourier components associated with signals which were ap- 
proximately constant over the duration of each data record, a short MATLAB routine 
was written to perform the fast Fourier transform of the signal and renormalize the 
large Fourier components to a given maximum intensity. Fig. 2 shows the fast Fourier 
transform of a typical RF signal before and after this procedure was applied. In each 
case the data were acquired using the "wide-band" filter configuration, whose response 
cuts off sharply below 23 MHz. 

The effect of digital filtering on detectability of a transient is illustrated in Fig. 3. The 
top panel shows the RF record whose Fourier transform was given in Fig. 2, on which has 
been superposed a simulated transient of peak amplitude 14.5 digitization units. (The 
data acquisition scale ranges from —128 to +127 digitization units; one scale division 
on the oscilloscope corresponds to 25 units.) The transient is invisible beneath the large 
amplitude associated with television and FM radio signals. The middle panel shows the 
result after application of the Fourier coefficient shrinkage algorithm. 

The event in Figs. 2 and 3 consisted of 32,768 data points obtained at a 1 ns sampling 
interval, with the trigger at the 20,000th point. The frequency resolution in the fast 
Fourier transform is thus 500 MHz (the Nyquist frequency) divided by 16,384, or about 
30 kHz. This permits rather fine distinction between frequencies containing a strong 
carrier and those which correspond to its weaker sidebands. At the same time, it permits 
time resolution to be preserved, allowing for the examination of rather rapid transients. 
To the extent that these transients do not contain Fourier coefficients exceeding a pre- 
determined threshold, they should be relatively unaffected by the shrinkage algorithm 
in the absence of interfering signals. However, since at Dugway signals in nearly the 
whole FM band (88-108 MHz) exceed the threshold, some distortion is unavoidable 
using such a method. In obtaining bounds on pulse amplitudes we therefore employ 
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Figure 2: Top panel: Fourier spectrum (in arbitrary units) of RF signals acquired at 
Dugway site using high-pass 25 MHz and low-pass 250 MHz filters. Prominent features 
include video and audio carriers for TV Channels 2, 4, 5, 7, and 11 (see Table 1 for 
frequencies), and the FM broadcast band between 88 and 108 MHz. Bottom panel: 
Fourier spectrum (same vertical scale) after renormalization of large Fourier components 
to a magnitude chosen here to be 3.16 x 10 3 . In practice best sensitivity to transients 
was obtained by renormalizing to a magnitude of 10 3 . 
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Figure 3: Effect of Fourier coefficient shrinkage on detectability of a transient. Top 
panel: raw RF record (in arbitrary units) with simulated signal superposed. Middle 
panel: record (same scale) after Fourier coefficient shrinkage. Here a maximum Fourier 
coefficient magnitude of 10 3 (in the units of Fig. 2) has been imposed. Bottom panel: 
the same record after denoising with a 10-point symmlet level L = 4 routine 13*o]. 
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a method involving the comparison of Fourier power in a given time window with the 
average power obtained over the whole data record for each Fourier component. This 
method is described below. 

4.5.2 Time-frequency analyses 

One can perform a fast Fourier transform using a small time window (typically 1024 ns) 
which is advanced sequentially through the data record, typically in steps of 100 ns. The 
frequency resolution of any given "snapshot" is then 500 MHz divided by (typically) 512, 
or a bit better than 1 MHz. A two-dimensional display of time vs. frequency then allows 
one to distinguish short transients (with components over many frequency bins) from 
continuous RF sources (with components in narrow frequency ranges over the entire time 
record). One such plot appears in Fig. 11, Sec. 5.3.1, below. In practice one may wish 
to suppress frequencies corresponding to the whole FM band and known TV stations, 
so as not to overload the dynamic range of the display. An alternative method [33] is to 
renormalize each point in time-frequency space so that deviations from the average in 
each frequency bin are displayed. This method is described further in Sec. 5.2. It was 
used for the main part of data analysis. 

4.5.3 Wavelet techniques 

The wavelet package Wavelab [34J contains a denoising routine which was adapted for 
our purposes. While an exhaustive search for optimized methods was not performed, 
good results in reducing noise levels were obtained using a 10-point symmlet routine 
with level L = 4 [33]. An example of a denoised signal is shown in the bottom panel 
of Fig. 3. Here a simulated signal of positive peak amplitude 14.5 digitization units has 
been added to an RF record otherwise free of transients. The effect of wavelet denoising 
is to reduce the amplitude of random high-frequency fluctuations while preserving edge 
effects such as transients. 

4.6 Signal simulation 

We wished to quantify the improvement associated with each method of signal process- 
ing. We thus simulated the expected signal by generating it using an arbitrary waveform 
generator, feeding it through the same preamplifier and filter configurations used for data 
acquisition, and superposing it on records otherwise free of transients. We successively 
reduced the amplitude of the superposed test signal until it could not be distinguished 
from random noise peaks, thereby obtaining an estimate of sensitivity. 

A Hewlett-Packard Arbitrary Waveform Generator was used to generate signals 
whose characteristics are illustrated in Fig. 4. These signals were taken to have the 
form f{t) = 6{t)At 2 (e~ Bt — Ce~ Dt ) with the coefficient C chosen so that f(i) has no 
DC component, and D corresponding to a long duration of the negative-amplitude com- 
ponent. For all pulses we chose D = B/20, so that C = (8000) -1 cancels the DC 
component. The Fourier components of the test pulse fall off smoothly with frequency. 
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Figure 4: Analytic depiction of typical pulse presented to filter-preamplifier configura- 
tion. Top panel: time dependence of pulse f(t) = 9(t)t 2 [e~ 0At — e~ a02 */8000](t in ns); 
bottom panel: Fourier spectrum of pulse (calculated analytically). In the top panel, the 
short bar above the pulse denotes S, the time difference between onset and maximum, 
while the longer bar below the pulse denotes A, the duration of the positive component. 

The initial t 2 behavior was chosen so that both the test pulse and its first derivative 
vanish at t = 0, as might be expected for a pulse from a developing shower. 

The simulated pulses are summarized in Table 6. Instead of quoting the value of A, 
we quote the maximum positive value of the pulse, both before and after filtration and 
preamplification. These values of Vpk reflect choices for convenience in display on the 
oscilloscope, and are otherwise arbitrary. 

The shape of the pulse of Fig. 4 is affected by preamplification and filtration as shown 
in Figs. 5 (broad-band) and 6 (narrow-band). The noise in these figures and the sharp 
feature at 125 MHz in Fig. 5 are associated with the system used to generate the test 
pulse, and the fact that the Fourier transform is taken over a much longer time than the 
duration of the pulse. 

Systematic studies of signal-to-noise ratios have been performed so far only for the 
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Table 6: Parameters of test signals. 5 is the time between pulse onset and maximum, 
while A is the duration of the positive component of the pulse. Vpk is the peak (positive) 
input voltage to the filter-preamplifier configuration. The letter after the peak voltage 
denotes (a) narrow-band (25-35 MHz) or (b) broad-band (> 25 MHz) configuration 
(see Sec. 14.4. 3J1 . Kut is the peak-to-peak amplitude of the pulse emerging from the 
filter-preamplifier configuration. S is the scale factor with which data were recorded on 
oscilloscope. 
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simulated pulses with S = 5 ns applied to a broad-band front end [(b) in Table 6]. The 
value of S is a measure of the distance R of closest approach of the shower core |3j ■ This 
choice corresponds to a typical distance R ~ 200 m. A typical pulse of this type gave a 
front end output of 21 mV peak-to-peak, acquired at an oscilloscope sensitivity of 5 mV 
per division. Each division corresponds to 25 digitization units, so the peak-to-peak 
range is about 104 digitization units, or slightly less than half the dynamic range (255 
units, or 8 bits). Positive and negative peaks are thus about 52 digitization units each. 

The stored test signal is then multiplied by a scale factor and added algebraically 
to a collection of RF records in which, in general, randomly occurring transients will 
be present. One then inspects these records to see if the transient can be distinguished 
from random noise. 

For the broad-band data we estimated that pulses with input voltages corresponding 
to about 1/5 the original test pulse amplitude can be distinguished from average noise 
(not from noise spikes!). Since the original test pulse had a peak value of 1.3 mV, this 
corresponds to sensitivity to an antenna output of about Vpk — 260 /iV. The ability 
to detect such a pulse with an effective bandwidth of about 30 MHz corresponds to a 
threshold sensitivity at the level of order 3 /iV/m/MHz (see Appendix B). 

Preliminary studies of simulated pulses applied to the narrow-band front end suggest 
a considerably poorer achievable signal-to-noise ratio, despite the expectation that the 
signal should have a large portion of its energy between 23 and 37 MHz. It appears 
difficult to detect a pulse from the antenna below about 0.7 mV, which for a bandwidth 
of 14 MHz corresponds to a threshold sensitivity of 7 /xV/m/MHz, not sufficient for our 
purposes. Studies of possible improvements of the analysis algorithm for the narrow- 
band data are continuing. 
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Figure 5: Test pulse of Fig. 4 after broad-band filtration (> 25 MHz) and preamplifica- 
tion. Top panel: time dependence of pulse; bottom panel: Fourier spectrum of recorded 
pulse for -20 fis < t < 12.768 pis. 
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Figure 6: Test pulse of Fig. 4 after narrow-band nitration (25-35 MHz) and preamplifi- 
cation. Top panel: time dependence; bottom panel: Fourier spectrum of recorded pulse 
for -20 /is < t < 12.768 fis. 
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Table 7: Triggers associated with CASA operation taken under various conditions. 



Front end 


Macintosh 


Dell 


Total 


Narrow-band 


5849 (139.78 h) 


1952 (48.97 h) 


7801 (188.75 h) 


Broad-band 


9603 (272.57 h) 


5416 (121.62 h) 


15019 (394.18 h) 


Low-frequency 





505 (17.67 h) 


505 (17.67 h) 


Total 


15452 (412.35 h) 


7873 (188.25 h) 


23325 (600.6 h) 



Table 8: Broad-band data recorded on Macintosh Quadra. 



Antenna 


CASA 


CASA 


Partial 


Total 


Polarization 


HV on 


HV off 


CASA HV 


events 


East-West 


4966 (119.03 h) 


859 (21.88 h) 


1957 (53.08 h) 


7782 (194.0 h) 


North-South 


696 (30.53 h) 


582 (23.25 h) 


543 (24.78 h) 


1821 (78.57 h) 


Total events 


5662 (149.57 h) 


1441 (45.14 h) 


2500 (77.87 h) 


9603 (272.57 h) 



5 Results 

5.1 Event sample 

More than 20000 triggers, obtained under various conditions of filtering, preamplifica- 
tion, and noise reduction during the period February 1997- March 1998, are summarized 
in Table 7. Events recorded on a Macintosh Quadra 950 and those recorded on a Dell 
LM Latitude laptop computer are listed separately because the power supply of the lat- 
ter introduced spurious transients. Our initial analysis concentrated on data taken with 
the Macintosh. For reasons mentioned above, we consider only the broad-band data at 
this time. Thus, our usable sample consists of over 9000 CASA triggers. In addition, 
periodic forced triggers were taken to monitor noise activity not associated with CASA 
operation. 

The broad-band data recorded on the Macintosh Quadra, summarized in Table 8, 
are subdivided into several categories. Data were taken with both East- West (EW) and 
North-South (NS) antenna polarizations. Moreover, since noise from CASA boxes was 
found to be a significant source of RF transients, data were taken with some or all CASA 
boxes disabled by turning off the high voltage (HV) supply to the photomultipliers. Even 
when HV is supplied only to boxes that are further than 100 m from the antenna, the RF 
transients from these boxes provide a strong background. (See the discussion in Sec. 15.21 
and Fig. 10 below.) One is unlikely to distinguish the RF pulses of the showers from 
this noise. Therefore, we concentrated on data with CASA HV off, with 859 triggers 
taken with EW antenna polarization (21.88 active hours) and 582 triggers taken with 
NS antenna polarization (23.25 active hours). For subsequent sensitivity calculations, a 
subset of data was used consisting of 756 EW triggers (17.25 hours) and 528 NS triggers 
(21.12 hours). The remaining data with CASA HV off occurred in very short runs (52 
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Figure 7: Top panel: intensity-vs.-time plot for maxima of 880 pulses recorded in 
698 triggers in January 1998 with CASA HV supplied to all stations. Bottom panel: 
time distribution of transients. All events recorded with East- West antenna polarization. 



EW triggers and 19 NS triggers) or was contaminated by local VHF communication 
signals (51 EW triggers and 35 NS triggers). Data with CASA HV on or partially 
disabled were not used for the present analysis. 



5.2 Characterization of transients associated with CASA op- 
eration 

Several means were used to characterize transients. One method with good time resolu- 
tion involved the shrinkage of large Fourier coefficients to a fixed maximum intensity, as 
in Figs. 2 and 3. Another, which we have used for results to be presented below, involves 
generation of a time- vs. -frequency intensity plot by Fourier-transforming 1024-ns subsets 
of the 50 /is data record, spaced by 100 ns steps. Since the data are sampled at 1 ns 
intervals, the frequency resolution of this method is thus about 1 MHz. The intensity 
S(u, t) is then averaged over time t for each frequency v to form an average intensity 
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Intensity vs. time, 1/98 transients, CASA HV off (824 entries) 
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Figure 8: Top panel: intensity-vs.-time plot for maxima of 824 pulses in 691 triggers 
recorded in January 1998 with CASA HV disabled. Bottom panel: time distribution of 
transients. All events recorded with East- West antenna polarization. 



S{v). The quantity S{v,t)/S{v) is an estimate of the degree to which the intensity at 
a given frequency v and time t exceeds the average over the 50 /is sampling time. We 
then average S(is, t)/S(v) over v to search for events in which the average intensity at a 
given time is exceeded in many simultaneous frequency bands. 

One can then search for peaks of each data record (there may be several peaks in a 
record), plotting intensity of their maxima against time relative to the trigger. One such 
plot is shown in Fig. 7 for a data run in which CASA HV was delivered to all boxes. A 
strong accumulation of transients, mostly with intensity just above the arbitrarily chosen 
threshold (mean + 3 a for each trigger sample), is visible at times —5 to —7 fis relative to 
the trigger. In a comparable plot for a run in which CASA HV was completely disabled 
(Fig. 8), only a small accumulation at times —6 to —7 /is is present. This excess appears 
due to transients with predominantly high-frequency components (over 100 MHz). Since 
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Figure 9: Signal of a typical transient associated with CASA operation. Top panel: 
before denoising; bottom panel: after denoising. 



signal pulses are expected to have more power below 100 MHz (see Fig. 5, bottom) we 
believe that this accumulation is not due to shower radiation, but most likely arises from 
the muon patches, one of which is within 75 m of the antenna. 

A typical transient occurring in a run with CASA HV on is shown in Fig. 9. The 
transients are highly suppressed (though not in all runs) when CASA boxes within 100 m 
of the antenna are disabled, as shown in Fig. 10. 

The time distribution of pulse maxima above an arbitrary threshold for 880 pulses 
detected with CASA HV on (one run from January 1998 composed of 698 files of data) 
is shown in the bottom panel of Fig. 7. The mean arrival time is about 6 /as before the 
trigger, with a distribution which is slightly broader for pulses arriving earlier than the 
mean. This broadening may correspond to some jitter in forming the trigger pulse from 
the sum of muon patch pulses. 

As mentioned earlier, the time for the trigger pulse to propagate from the central 
station to the RF trailer was measured to be 2.15 lis. One expects a similar or slightly 
greater travel time for pulses to arrive from muon patches to the central station (see 
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Intensity vs. time, 1/98, CASA HV on > 100 rn from ant. (903 entries) 
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Figure 10: Top panel: intensity-vs.-time plot for maxima of 903 pulses recorded in 620 
triggers in January 1998 with CASA HV disabled for boxes within 100 m of antenna. 
Bottom panel: time distribution of transients. All events recorded with East- West 
antenna polarization. 

Fig. 1). Moreover, the muon patch signals are subjected to delays so that they all arrive 
at the central station at the same time for a vertically incident shower. Thus, the peak 
in Fig. 7 is consistent with being associated with the initial detection of a shower by 
CASA boxes. This circumstance was checked by recording CASA trigger request signals 
simultaneously with other data; they coincide with transients such as those illustrated 
in Fig. 9 within better than 1/2 /is. 

The RF signals from the shower are expected to arrive no later than, or at most 
several hundred nanoseconds before, the transients associated with CASA operation. 
They would propagate directly from the shower to the antenna, whereas transients from 
CASA stations are associated with a longer total path length from the shower via the 
CASA station to the antenna. There will also be some small delay at a CASA station 
in forming the trigger request pulse. Thus, we expect a genuine signal also to show up 
around 6-7 /is before the trigger. 

The time coincidence of the CASA RF transients and the shower signals is a signifi- 
cant obstacle to detecting genuine pulses. Therefore, data with CASA HV on or partially 
disabled were not used for the present analysis. As we show below in Sees. 5.3.1-5.3.3, 
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no significant peak is visible around 6-7 fis before the trigger for data recorded with 
CASA HV off. The upper limit on the rate of events giving rise to such a peak can 
be used to set a limit on RF pulses associated with air showers, as we demonstrate in 
Sec. 5.3.4. 

5.3 Estimated upper bounds on broad-band signals 

As mentioned in Sec. 5.1, the following discussion is based on 17.25 active hours of 
data accumulation with EW antenna polarization and 21.12 hours with NS antenna 
polarization. The small duration of this subset of data limits its sensitivity to RF 
signals from the shower. 

5.3.1 Criteria used to distinguish noise and signal transients 

The main difficulty associated with pulse detection is that signal pulses are not easily 
distinguishable from large spurious pulses originating from atmospheric discharges. Both 
air shower pulses and these background noise pulses can considerably exceed the average 
noise level. Several criteria can be used to distinguish signal pulses from noise. The 
conventional criteria of the previous studies have been that the pulse should be (1) 
larger than the average noise level by some small specified amount, (2) time coincident 
with shower particles, and (3) bandwidth limited All these criteria were adopted in 
this study and one more has been added: The pulse should have approximately uniform 
distribution over frequency within its limited bandwidth [see (c) below]. 

We now note the particular criteria used to distinguish signal pulses in this study. 

(a) Pulse magnitude 

Continuous RF interference in each frequency channel was removed by advancing a 
moving 1024-ns window in 100 ns steps through the data record to produce a time-vs.- 
frequency plot [Sec. I4.5.2] and then using the averaging procedure described in Sec. 5.2. 
Defining the average signal as 1 (in arbitrary units), the pulse threshold in each event 
was taken to be the larger of either (a) the mean plus three standard deviations, or (b) 
1.8 (in the same units). The former permitted removal of an average noise level; the 
latter discriminated against small noise transients. The final result was not affected by 
the choice of the factor 1.8 since the subsequent analysis [Sec. I5.3."3~] used a considerably 
higher threshold. 

(b) Limited bandwidth of pulses 

Broad-band filtering limits the frequency range to 23-250 MHz [Sec. I4.4.3j . The 
investigation of the simulated pulses on a 2d-plot of intensity vs. time and frequency 
suggested that, unlike some strong noise transients, the signal pulse intensity declines 
drastically in the range above approximately 100 MHz (Fig. 11). This feature is consis- 
tent with theoretical predictions for the shower pulse spectrum [HE]. R can be chosen 
as a criterion for separating noise and signal transients. The ratio of mean intensities 
averaged over 23-100 MHz relative to that averaged over 100-250 MHz was found to be 
greater than 1 for all simulated pulses and smaller than 1 for some noise pulses. 

One can consider intensities averaged over the part of the whole frequency range. 
To facilitate separation of signal transients it is preferable to choose a region where 
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Figure 11: Simulated pulse (horizontal band at fis) and noise transient (horizontal 
band at about 15 us) on a 2-d plot of intensity (grayscale) vs. time (vertical axis) and 
frequency (horizontal axis). Unlike the simulated pulse, the noise transient contains high 
frequency components. Grayscale: black color denotes the highest intensities, white the 
lowest. Vertical bands indicate continuous RF sources in the 54-200 MHz range. 



the signal-to-noise ratio is particularly large. Unfortunately, the whole region from 
23 MHz to 100 MHz cannot be effectively used for this purpose. The 54-82 MHz range 
is occupied by TV channels 2, 4 and 5, which leads to high noise levels. The same is 
true for the whole FM band (88-108 MHz) [Sec. I4.2.1j . However, this is not the case in 
the 24-54 MHz range. The noise level in this range is mostly uniform, and simulated 
pulse intensities are particularly large there in comparison with the noise level. The 
assumption that the 24-54 MHz range provides the best signal-to-noise intensity ratio 
has been tested. The results for this range have been compared with the ones obtained 
in the 10-54 MHz and 24-86 MHz regions and were found to be superior. Subsequently, 
the range of 24-54 MHz was chosen as the main region of investigation. 

After that, it was natural to choose the ratio of mean intensities averaged over 24- 
54 MHz relative to that averaged over 55-250 MHz at the moment of each pulse as a 
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Table 9: Fraction of noise pulses passing criterion (a) that can also pass criteria (c) and 
(b) with ratio parameter threshold of 1.4, and whose maximum intensities are larger than 
some very high intensity threshold. Fraction of simulated pulses of different strengths 
(in /iV/m/MHz), that can pass criteria (c) and (b) with ratio parameter threshold of 1.4, 
and whose maximum intensities (after they are superimposed on noise) are larger than 
the same high intensity threshold. The strengths of the pulses are chosen to be stronger 
than 1/n times an original test pulse amplitude of 13.7 //V/m/MHz, n — 2,3,4, ... ,7. 





Noise pulses 


Simulated pulses that are stronger than 


6.85 


4.57 3.43 2.74 


2.28 


1.96 


/iV/m/MHz 


North-South 


4.47% 


94.1% 


91.7% 


85.2% 


66.9% 


47.7% 


36.7% 


East- West 


4.75% 


98.8% 


93.4% 


84.1% 


67.7% 


48.7% 


41.4% 



criterion for discriminating noise and signal pulses. For all simulated pulses this ratio 
was greater than 1. All pulses for which this ratio was smaller than 1 were assumed to 
be noise transients and discarded. Also, a ratio parameter threshold other than 1 can 
be chosen. The parameters that provided best results were found to lie between 1.4 and 
1.8 [Section EH - 

(c) Approximately uniform distribution of pulse intensity over frequency in the 24~ 
54 MHz range 

The most intense noise transients that met criteria (a) and (b) were found to display 
a peculiar feature: Their intensities were concentrated in a small region of approximately 
10 MHz width somewhere in the 24-54 MHz range. This non-uniformity allowed such 
pulses to be ruled out. If the pulse intensity, integrated over any 9 MHz width window 
(in the 24-54 MHz range), was greater than the intensity integrated over the remaining 
21 MHz, then such a non-uniform pulse was discarded as a noise transient. The reason 
for choosing a 9 MHz width window was that most simulated pulses passed this test, 
while many noise transients did not. 

The effectiveness of these three criteria is illustrated in Table EJ For simulated 
pulses that are stronger than 3.43 //V/m/MHz the combination of criteria (b), (c) and a 
very high intensity threshold provides a 84.1/4.75 ~ 18 times increase for their relative 
fraction with respect to noise transients. 

5.3.2 Method outline 

The time coincidence of the pulse with air shower particles implies that a signal pulse 
should be looked for in the range between —7 and —6 fis [Sec. 15. 2j . which will be 
referred to in what follows as "the interesting time bin" . It is important to note that the 
aforementioned test criteria are not efficient if applied only to the pulses found in this 
time bin. Indeed, some noise pulses in the bin meet all criteria. Hence, accepting all 
the pulses that pass this test would not guarantee that signal pulses are present among 
them at all. We developed a different approach. 
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Transients that meet criterion (a) 




°30 -20 -10 10 20 30 

Time of maximum relative to trigger (microseconds) 

Figure 12: Intensity vs. time plot for maxima of transients in 756 files of data with East- 
West antenna polarization and CASA HV off. The horizontal line denotes the threshold 
taken to lie at the level of the 17th strongest pulse. Only criterion (a) of Sec. 5.3.1 was 
employed, resulting in a large (1367) number of entries. 

The time distribution of noise pulses is assumed to be uniform. Hence, the following 
technique can be adopted. First, the time distribution histogram of the pulse maxima 
(from the whole data run) can be plotted. Then, in the absence of signal pulses, the 
number of pulses in a time bin should obey a Poisson distribution with the average being 
equal to the average pulse number per bin. Suppose that an accumulation of pulses in 
the interesting time bin (—7, —6) /is is large enough that its probability according to 
the Poisson distribution is less than 5%. Then this signifies the presence of signal pulses 
with a confidence level of 95%. 

Thus, a sufficiently large relative accumulation of entries in the interesting time bin 
was adopted as a key criterion in search of signal pulses. 
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5.3.3 Detailed description 

A routine was designed to scan all events in the whole data run and select the transients 
that passed criterion (a) in Sec. 5.3.1. The ratio parameter (see above) and uniformity 
parameter [1, if criterion (c) was satisfied, and 0, if not] were recorded along with time 
relative to the trigger and intensity of a transient. The intensity vs. time plot of pulse 
maxima selected in this manner for the data taken with East- West antenna polarization 
is shown in Fig. 12. 

In order to be detectable, signal pulses should be stronger than the average noise 
transient, whose intensity equals 2.5 as a result of the specific criteria imposed. Hence, 
it is reasonable to set the intensity threshold sufficiently high to enhance the relative 
accumulation in the (—7, —6) fis time bin. One can set a high threshold and determine 
the number of pulses that pass through it. Divided by number of bins, this number gives 
the average pulse number per bin. This, in turn, can be used as an average value of the 
Poisson distribution to determine whether the probability of the observed accumulation 
in the interesting bin is less than 5%. 

The effectiveness of this technique depends on the choice of intensity threshold. If the 
threshold is too low, the histogram includes many noise pulses. Then, the signal pulse 
accumulation in the interesting time bin becomes relatively small to be distinguished 
from statistical fluctuations. If the threshold is too high, it may cut out some signal 
pulses and the remaining ones may not make up a significant accumulation. The choice 
of the optimum threshold is discussed below. 

Optimum intensity threshold. Suppose we look for an accumulation of 2 or more 
entries. Such an accumulation is significant (i.e. its probability is < 5%) if the average 
number of pulses per bin is about 0.355. Indeed, according to Poisson distribution, 
P(0 entries) = exp(-0.355) « 0.701, P(l entry) = 0.355 • exp(-0.355) « 0.249, so 
P(> 2 entries) = 1 — 0.701 — 0.249 = 0.05. Since the standard histogram used in the 
study contained 48 bins (of 1 /is each), the average value of 0.355 pulses per bin would 
correspond to a total number of 0.355 • 48 ^ 17 pulses. These considerations suggest 
that the intensity threshold would not have some fixed value but would correspond to 
the intensity of the 17th strongest pulse [see Fig. 12]. Then, the total number of entries 
in the histogram is 17 and the average number of entries per bin is about 0.355. So, 2 
or more entries in the interesting bin, if observed, would have the probability of 5% and 
would reveal the presence of signal pulses with 95% CL. 

Naturally, the probability of 3 or more entries would be even smaller and approxi- 
mately equal to P(> 3 entries) = P(> 2 entries)-P(2 entries) = 0.05-^^ exp(-0. 355) 
~ 0.006. That is, if such an accumulation were observed, the presence of signal pulses 
could be claimed with a confidence level of 99.4%. However, in reality it is unlikely that 
signal pulses could be detected with such a confidence level. To improve the detection 
chances, it is advantageous to lower the threshold to the 39th strongest pulse. This sets 
the average number of entries per bin to 39/48 ~ 0.813 and makes the probability of 
3 entries equal to 5%. The threshold of the 39th strongest pulse is the lowest one that 
still guarantees a confidence level of at least 95%. 

Thus, the convenient threshold depends on the significant accumulation number one 
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Table 10: Significant accumulation of entries in the interesting time bin and the corre- 
sponding optimum threshold (calculated according to Poisson distribution). 
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is looking for. For any particular accumulation, one can choose the optimum threshold 
which would guarantee a confidence level of at least 95%. The significant accumulation 
number and the corresponding optimum threshold can be found in Table 

This method did not reveal any significant accumulation in the interesting time 
bin. So far, however, only criterion (a) and intensity thresholds were applied to de- 
tected pulses. As was already mentioned in the beginning of this Section, the ratio and 
uniformity parameters were recorded for each simulated or noise pulse. This made it 
easy to impose criterion (b) with different ratio parameter thresholds and criterion (c). 
Their application should have increased the ratio of signal pulse number to noise pulse 
number. Unfortunately, no significant accumulation has been detected (see Fig. 13). 
This means that signal pulses are quite rare, so that the resulting accumulation is not 
significant. However, upper limits can be placed on the rate (in h -1 ) of signal pulses 
that are stronger than some fixed value. This will be our aim for the next subsection. 
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Time distribution histogram of the pulse maxima (65 entries) 

4 1 i i i i ■ i i i i 



c/) 
CD 



CD 



CD 
-O 

E 

13 



■ ■ ■■■ ■ ■■■■■ ■ ■ ■ i ■ ■■■■■■■■■■■■■ ■■ ■ ■■■■■■ 

—20 -15 -10 -5 5 10 15 20 

Time relative to trigger (|j,sec) 

Figure 13: Time distribution histogram for maxima of the 65 strongest pulses extracted 
from 756 files of the EW data with (a), (b) and (c) imposed. The ratio parameter 
threshold for criterion (b) is 1.4. An accumulation of 4 or more pulses in the (-7,-6) /j,s 
bin would signify the presence of signal pulses. Such an accumulation was not detected. 

5.3.4 Upper limits on the rate 

The total number n of observed transients in the interesting time bin consists of both 
noise and signal events. The former are Poisson-distributed with known average /i. The 
latter obey a binomial distribution with a known probability P to pass thresholds and 
an unknown total number iV of signal events contained in the run. Using a unified 
approach to the statistical analysis of small signals [HZ], we construct 95% confidence 
belts for unknown N. For any choice of imposed thresholds, the lower end of confidence 
intervals is 0, i.e. only upper limit on the total number of signals in the run can be 
placed. 

Probability P s (k) that the interesting bin contains k signal events. For- 
mula (1) [Section 3] shows that £ u is directly proportional to primary energy E p . The 
differential rate of primaries is known to fall approximately as 1/-E- [21] • Similarly, the 
differential rate of £„s should be proportional to 
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Simulated pulses of different strength £ v greater than some £® can be added to each 
event of the data run. To simulate the expected rate for £ u , a Monte-Carlo simulation 
was performed in such a way that the number of added simulated pulses with strength 
£ v falls as l/£%. 

Each event underwent the averaging procedure described in Sec. 5.2. The averaged 
intensities of simulated pulses were compared to the values of the 17th (or 39th, 65th, 
etc.) strongest noise pulse of the initial data (without added simulated pulses). The 
fraction of simulated signal pulses that were higher than this threshold gives the proba- 
bility P for a signal pulse with strength greater than some £® to exceed this threshold. 
We determined the value of P as a function of £®, the intensity threshold, and the set 
of criteria ((a), (b),(c), or some combination of them) imposed on both noise and signal 
pulses. Then, these values were used to make an analytical estimate of the probability 
to detect a significant accumulation number of entries in the interesting bin. 

By making an assumption for the total number of signal pulses N during the run and 
using the binomial distribution, one can calculate the probability that k out of these N 
signal pulses exceed the threshold: P s (k) = C(N, k) P k (1 — P) N ~ k , where C(N, k) is the 
standard binomial coefficient. This probability not only depends on £®, the intensity 
threshold, and the set of imposed criteria but also on the aforementioned assumption 
for the number N. 

Determining upper limits N up . Suppose, for instance, that we are looking for 3 
or more significant entries in a bin. 

First, we impose a particular set of criteria on all pulses in the run, pick the 39 
strongest of them and determine the number n of entries in the interesting bin. 39 
pulses in the histogram correspond to an average /z = 39/48 ~ 0.813 entries per bin. 
The Poisson distribution P n {k) = /i fc e~ M /fc! determines the probability to have k pulses 
in a bin. 

Second, we set the intensity threshold at the level of the 39th pulse. Then we 
determine P s (k) for the same set of imposed criteria, this intensity threshold and any 
total number iV of signal transients in the run. 

Consider the construction of an acceptance interval of n values for some fixed total 
number N of signal transients in the run. n entries in the interesting bin can be the 
result of different combinations of noise and signal entries. We calculate the probability 
to detect n entries in the bin as 

n n 

P(n\N) = P s (k)P n (n - k) = e - " £ C{N, k) P k (1 - P) N ~ k ii n - k /(n - k)\ 

k=0 k=0 

Take some values of n and N, for example, n = no and N = Nq. The probability 
P(n \N ) might be small but not so small with respect to P(rio\N best ), where N best is such 
an alternate hypothesis which maximizes P(uq\N). The ratio R = P(nQ\N Q )/P(nQ\Nb est ) 
is the basis of the ordering principle outlined in [3~7| . 

For any N we add values of n into the acceptance interval in decreasing order of 
R. As soon as the sum of P(n\N) exceeds the confidence level of 0.95, the acceptance 
interval is completed. The confidence belts that can be constructed using this procedure 
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Figure 14: Confidence belt based on ordering principle of reference [37], for 95% CL 
intervals for unknown total number N of signal pulses in the run. The probability 
P = 0.841 for a signal pulse to go through a set of cuts, and the Poisson background 
mean fi = 65/48 = 1.354, are the parameters of the plot. The presence of n > 4 entries 
in the interesting time bin from the total of 65 entries in the histogram would signify the 
presence of signal pulses in the run [see Table 10]. This is also reflected in this Figure 
as for n > 4 N low becomes nonzero. P = 0.841 corresponds to the probability for signal 
pulses stronger than 1/4 the original test pulse amplitude to pass criteria (b), (c), and an 
intensity threshold at the level of the 65th strongest pulse of the EW data [see Table 9]. 
The ratio parameter threshold for criterion (b) is 1.4. Only one of the 65 strongest pulses 
was observed in the interesting time bin [see Fig. 13], implying the upper limit N up = 2 
(see also the corresponding dot in Fig. 15(b) at £® = 3.43 /iV/m/MHz). 



are shown in Fig. 14. For any observed number n of entries in the interesting bin, these 
belts provide intervals for allowed values of N. For any choice of imposed criteria and 
for all values of n measured in the experiment, the lower end of the intervals was 0. 
Thus, the actual signal transients were not detected and only the upper limit N up on 
their number during the run can be placed. 

Upper limit R up {£ v > £®). Divided by the total run time (17.245 hours with EW 
antenna polarization and 21.116 hours with NS polarization), N up gives the upper limit 
on the rate of signal pulses of strength greater than £°. We will denote it R up {£ v > £®). 
Of course, one can obtain different values of R up (£ u > £®) when looking for a different 
significant number of entries in the interesting bin. We searched for the accumulation 
values from 2 to 20 in the bin (—7, —6) fis. The lowest R up (£ u > £®) gives the most 
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stringent upper limit on the rate. Its value also depends on criteria imposed on both 
signal and noise pulses. The best results were achieved when all three criteria (a), (b) 
and (c) were employed. The ratio threshold parameter for criterion (b) was tested in the 
region from 1 to 2.4 with step 0.2. The thresholds that give the best results were found to 
be 1.4 for EW and both 1.6 and 1.8 for NS data. Compared to the analysis with neither 
(b) nor (c) employed, these threshold parameters provide up to 28% lower upper limits. 
The resulting values of R up (£ u > £°) are shown as dots in Fig. 15 for several values of 
£® for the data taken with both East- West and North-South polarizations. Signal pulses 
stronger than 4.57 /xV/m/MHz are expected to be quite rare and were definitely absent 
from the NS data, just like pulses stronger than 6.85 //V/m/MHz were absent from the 
EW data. Zero upper limits on the rates of these events are indicators of these facts. 

As was mentioned above, the differential rate R{£ v ) should be proportional to l/£%. 
Then, the rate of signal pulses of the strength greater than £® is 



i.e., proportional to l/(£°) 2 . Hence, the lowest curve cj (£°) 2 passing through one of the 
nonzero dots (solid line in Fig. 15) represents the strictest upper limit on this rate. The 
value of the coefficient c was found to be 0.555 for the NS data and 0.889 for the EW 
data: 



where £® is in //V/m/MHz. These results should be interpreted as the upper limits 
on the rates of events where the North-South (or East- West) projection of the signal 
pulse was greater than some value. We will use them in Appendix A for sensitivity 
calculations. 

5.4 Discussion and summary 

The small duration of the present subset of data prevents us from confirming or refuting 
previous claims for shower pulses. As mentioned, this subset was taken with CASA HV 
off, and was initially collected in order to check indications of a signal which occurred 
with CASA HV on. Using the estimate of 2.25 strong signal pulses per day [Sec. 4.4.4], 
one expects about three or four detectable shower pulses during 38.36 hours of obser- 
vation. Three or four detectable particles that might give a strong signal during the 
experiment would be too scanty an amount to be detected over the noise background, 
which closely mimics genuine signals. 

Nonetheless, we were able to place upper limits on the integral rate of signal pulses. 
Under certain simplifying assumptions these upper limits can be used for evaluating the 
numerical calibrating factor s. This factor was set to 20 in Eq. but there remained 
an uncertainty regarding its value. In fact, subsequent reports claimed values as small 
as 1.6. We could only set upper limits on s: s < 31 from the EW and s < 34 from the 
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Strength E y ° ((iV/m/MHz) 

Figure 15: Upper limit R up (S u > £®) on the rate of signal pulses stronger than a fixed 
strength £®. Plotted dots show the upper limit for antenna system gain G — 1.8. Small 
numbers near the dots indicate N up , the upper limit on the number of the signal events 
in the run. The dots were plotted for the pulses that are stronger than 1/n times 
the original test pulse amplitude of 13.7 /iV/m/MHz, n = 2,3, . . . , 7. The solid line 
represents the lowest upper limit corresponding to them. Data were taken with (a) 
North-South (21.12 hours) and (b) East- West (17.25 hours) antenna polarization. 
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Figure 16: Rate of events at the CASA detection area in which the East-West projection 
of the signal pulse is greater than some £ °. The lower three curves show the predicted 
rates based on different values of s: s = 20 ± 9.5 (dash-dotted line, the original Haverah 
Park result [3] based on ~100 detected shower pulses); s = 9.2 ±1.4 (dashed line, 
the Soviet group result jlj based on ~1000 pulses), and s = 1.6 ± 0.24 (dotted line, 
the updated Haverah Park result |3J based on ~1000 pulses). The uncertainty in the 
latter two values is 15% jl]. We evaluated the uncertainty in the former as being y/l0 ■ 
15% ~ 47%. In this logarithmic plot the error bars have constant lengths along their 
corresponding curves. The top curve is the upper limit in this paper which corresponds 
to s = 31. These four curves are determined under the assumption of a transverse 
current radiation mechanism [TU] . 
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Figure 17: Rate of events at the CASA detection area in which the North-South pro- 
jection of the signal pulse is greater than some £®. The lower three curves show the 
predicted rates based on different values of s: s = 20 ±9.5 (dash-dotted line, the original 
Haverah Park result |3] based on ~100 detected shower pulses); s = 9.2 ± 1.4 (dashed 
line, the Soviet group result jl] based on ~1000 pulses), and s — 1.6 ± 0.24 (dotted 
line, the updated Haverah Park result jl] based on ~1000 pulses). The uncertainty in 
the latter two values is 15% j3j. We evaluated the uncertainty in the former as being 
VTO • 15% ~ 47%. In this logarithmic plot the error bars have constant lengths along 
their corresponding curves. The top curve is the upper limit in this paper which corre- 
sponds to s = 34. These four curves are determined under the assumption of a particular 
radiation mechanism discussed in Appendix A. 
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NS data [see Appendix A]. The former was obtained under the assumption of the trans- 
verse current radiation mechanism, the latter under an alternative radiation mechanism 
discussed in Appendix A. Figs. 16 and 17 compare our results with those of previous 
experiments. 

Had additional observation time with CASA boxes disabled been available, more 
signal pulses and more noise transients would be recorded. Would it be easier to prove 
the presence of shower signals in the data? Taking into account that the average signal 
pulse from a high energy air shower is much stronger than the average noise pulse, the 
content of the most intense transients (those stronger than 17th strongest, 39th, etc.) 
would shift in favor of more signal pulses. Thus, additional running time with CASA 
HV off would have allowed one to at least place a stricter upper limit on the rate of 
shower pulses, if not to detect them. 

6 Issues specific to a giant array 

6.1 Frequency range 

The study of pulse components above 30 MHz remains a useful restriction in view 
of the changing RF environment encountered at lower frequencies. During years of 
sunspot minima the whole range of frequencies above 23 MHz remained relatively clear 
of broadcast interference, while as the sunspot numbers increased toward the end of 
1997, daytime interference became particularly intense from the citizen's band just above 
27 MHz. This interference typically subsided at sunset. If observations below 30 MHz 
are contemplated, they would be most useful during years of expected sunspot minima 
(e.g., 2006-7), to minimize effects of long-distance ionospheric reflections. 

6.2 Number and spacing of receiving sites 

It is expected |3] that signals above 30 MHz decrease rapidly as a function of distance 
between the antenna and the shower core's closest approach. Thus, we expect that in 
the Auger array, with spacing of 1.5 km between sites in a hexagonal array, RF stations 
would have to be distributed with roughly the same or greater density. One possibility 
for minimizing RF interference from Auger stations would be to place the RF station 
at the interstices between them. This would raise the cost per station, since it would 
require microwave communication with Auger sites and auxiliary sources of power. 

At each station it may be helpful to have two antennas, one registering pulses of 
east-west polarization and one for north-south polarization. Differential signals as well 
as individual ones should be recorded. Coincidences among several stations may be 
associated with particularly large showers. Frequency-dependent antenna phase response 
should be modeled or measured, as noted in Sec. 4.4.8. 

6.3 Digitization requirements 

Previous work by one of us (J.F.W.) involved detection of electromagnetic pulses, in- 
cluding those possibly produced by cosmic-ray-induced electromagnetic discharges, with 
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frequencies in the 30 - 100 MHz range. Part of this work included building hardware 
for self-triggering on short duration wide-band RF pulses. Many of the pulse identifica- 
tion, fast- digitization and memory problems were identical to those for pulse detection 
at CASA/MIA. Time-frequency plots were obtained similar to those one would generate 
in a survey at CASA/MIA. Similar requirements also are encountered for digitization of 
data from the KamLAND Experiment [3*%] . 

Our experience with the present system indicates the need for expanded dynamic 
range if continuous-wave sources of RF interference (such as FM and TV broadcast 
stations) are to be eliminated digitally. Thus, one needs at least 10-bit and probably 12- 
bit range, with a digitization rate of at least 400 MHz so as to be sensitive to frequencies 
up to 200 MHz. Although the expected signal is likely to be concentrated at lower 
frequencies (probably below 100 MHz), the expanded frequency range has proved useful 
in distinguishing expected signals from other transients. 

6.4 Stand-alone trigger 

Our results do not indicate that a trigger based on RF signals alone can yield useful 
correlations with air shower events. This result may be specific to the location of the 
CASA/MIA array; such a trigger may be less subject to noise at a remote location 
such as the Southern Hemisphere Auger site. An RF survey performed there would 
be useful in determining the utility of such a trigger. At such a site, more free from 
man-made noise than the CASA/MIA site, one would have to perform further studies 
allowing discrimination between random triggers (such as those induced by atmospheric 
discharges) and those induced by air showers. One would also attempt to detect galactic 
noise as a further indication that the site was sufficiently quiet. 

6.5 Integration into the data stream 

The data of CASA/MIA and that of the RF detection experiment were only integrated 
off-line. Any further studies should allow for simultaneous acquisition of both sets 
of data. Since the Auger project proposes to use microwave communication between 
stations, this same link should be considered for communicating RF signal acquisition 
results to a central data stream. 

6.6 Status of GHz detection 

David Wilkinson, who visited the University of Chicago during the spring of 1995, has 
proposed looking into the power radiated at frequencies of several GHz, where new op- 
portunities exist associated with the availability of low-noise receivers. These techniques 
have now been implemented in the RICE project |39 4 , which seeks to detect pulses with 
frequency components around 250 MHz in Antarctic polar ice. 
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Table 11: Television stations within 400 km of Dugway on channels not assigned to the 
Salt Lake City metropolitan area. Ref. also lists four "new" (unidentified) stations 
on Channel 3 for Price, UT (distance ~ 160 km, heading ~ 110 degrees) and two on 
Channel 12 for Logan, UT (distance 178 km, heading 14 degrees). 



Call 


Location 


Channel 


Power 


Distance from 


Heading 


sign 






(kW) 


Dugway (km) 


(degrees) 


KBJN 


Ely, NV 


3 


100 


212 


239 


KIDK 


Idaho Falls, ID 


3 


100 


363 


1 


KBNY 


Ely, NV 


6 


100 


212 


239 


KPVI 


Pocatello, ID 


6 


100 


301 


6 


KBCJ 


Vernal, UT 


6 


83.2 


306 


87 


KIFI-TV 


Idaho Falls, ID 


8 


316 


363 


1 


KENV 


Elko, NV 


10 


3.09 


272 


280 


KISU-TV 


Pocatello, ID 


10 


123 


301 


6 


KUSG 


St. George, UT 


12 


10 


359 


191 



6.7 Other options 

Dispersion between arrival times of GPS signals on two different frequencies may serve 
as a useful monitor of air shower activity. The possibility of correlation of large showers 
with such dispersion events could be investigated. 

It may be possible to monitor commercial broadcast signals in the 54-216 MHz range 
to detect momentary enhancements associated with large showers, in the same sense that 
meteor showers produce such enhancements. Television channels for which no nearby 
stations exist offer one possibility. The data taken at CASA/MIA have not yet been 
analyzed in terms of such enhancements, but represent a potential source of information. 
In Table we note the locations of TV stations broadcasting on VHF channels other 
than those assigned to the Salt Lake City metropolitan area within 400 km of Dugway. 
These are channels 3 (60-66 MHz), 6 (82-88 MHz), 8 (180-186 MHz), 10 (192-198 MHz), 
and 12 (204-210 MHz). The availability of at least two stations on Channel 3 and three 
on Channel 6 at greatly differing headings indicates that this method may have some 
promise. 

Radar detection of showers offers another exciting possibility 19 . This method 
resembles the use of distant fixed VHF stations for generating reflections off showers, 
but allows for a more carefully controlled environment. 

7 Conclusions 

A prototype system for the detection of radio-frequency (RF) pulses associated with 
extensive air showers of cosmic rays was tested at the Chicago Air Shower Array and 
Michigan Muon Array (CASA/MIA) in Dugway, Utah. This system was under con- 
sideration for use in conjunction with the Pierre Auger Project, which seeks to study 
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showers with energies above 10 19 eV. 

The system utilized a trigger based on the coincidence of 7 out of 8 buried muon 
detectors around the periphery of the CASA array Transients were indeed detected 
in conjunction with large showers, but they were identified as arising from the CASA 
modules themselves, most likely from the electronics generating trigger request (TRQ) 
pulses. Such transients could be eliminated when the high voltage (HV) on CASA 
phototubes was turned off; in such cases the muon trigger continued to function. By 
comparing upper limits on detected transients with simulated pulses, it was possible to 
place upper bounds on the rate of detection of RF pulses of various intensities. These 
upper bounds are summarized in Fig. 15; they typically involve rates of one every few 
hours for the largest field strengths claimed in the literature @] . Based on our estimates, 
it is unlikely that the present experiment can reach the sensitivity limits of the Haverah 
Park results, which reported lower field strengths in their latest work [4J. 

A number of lessons have been learned from the present exercise. These are probably 
most relevant for any installation contemplated in conjunction with the proposed Auger 
project [T3j . 

(1) One must take special care to survey transients produced by components of the 
array. For the Auger detector, one must install one or more antennae close to the 
proposed Cerenkov detectors and their associated digitizers, and study the response to 
artificially induced signals. Based on our experience, in which the present bounds are 
based on a small subset (38 hours) of the total data sample (600 hours), one should 
focus as soon as possible on a configuration in which usable data can be gathered. 

(2) The method of communication between Auger modules will affect what form of 
RF detection is feasible. If radio links employing microwave (> 800 MHz) frequencies 
are used, the present system will not be as seriously compromised as it apparently was 
by the communication system used at CASA/MIA. On the other hand, detection of RF 
signals above 1 GHz will suffer interference from such a system. 

(3) Consideration should be given to placement of antennas at sites sufficiently far 
from surface detector modules that pulses from these modules do not constitute a serious 
source of interference. In the Auger case, the modules are arranged in a triangular array 
with 1.5 km spacing, so that the maximum spacing between an antenna inside the array 
and any module could be as large as 1.5/v3 ~ 0.87 km if the antenna is placed at an 
equal distance from the three closest modules. 

(4) Coincidence between RF signals detected at several antennas is desirable, as was 
found in the earliest experiments j2J. 

(5) The digital filtering algorithms employed in the present study, although not yet 
pushed to their optimal efficiencies, appear to be limited by the 8-bit dynamic range 
employed in detection using a digitizing oscilloscope. Consideration should be given to 
a system with larger dynamic range, at least 10-bit but preferably 12-bit. 

(6) One can probably afford to economize by reducing the sampling rate, certainly 
to 500 MSa/s but perhaps as low as 200 MSa/s, since transients are expected to have 
their main frequency components below 100 MHz, and by reducing the active sampling 
window from the present value of 50 fis to a lower value, depending on the geometry of 
the array. 

(7) The frequency range studied in the present work (23-200 MHz) will be more 
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useful if continuous RF sources, such as FM and television stations, are much weaker 
than they were at the Dugway site. This is a possibility in the remote Argentine site at 
which the first Auger array is to be constructed [14] : a survey of field strengths there 
would be desirable. 

(8) One should take data simultaneously with two antennas polarized in perpendic- 
ular (EW and NS) directions. This allows the determination of two components of the 
electric vector, not just one of its projections. 

Although no RF signals have been detected in conjunction with CASA/MIA events, 
the present study has revealed a number of useful criteria for future experiments of 
similar type. It is hoped that a prototype at an Auger site will further focus these 
criteria. 
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Appendix A: Sensitivity Calculation 

One of the aims of this study was to evaluate the antenna response as a function of 
shower parameters. The dependence on R, 9, a, E p has been established in |3] (Eq. (1)). 
However, there remained an uncertainty as regards the calibrating factor s, for which 
experimental data do not supply an exact value. Indeed, we have quoted at least three 
versions of it: first, in Eq. (1) it is set to 20; second, the Haverah Park group subsequently 
reported that = 0.6, while the Soviet group claimed the value of 3.4 jlj. To make 
the connection between E^ and s, let us note that s ~ E^ / exp(— 1). Thus, we can 
infer that the Haverah Park group updated Eq. (1) to s — 1.6, and the Soviet group to 
s = 9.2. 

Ideally, we would like to calculate the calibrating factor precisely from the rate 
R(E u ew > £„)■> but unfortunately this rate cannot be obtained accurately from experi- 
ment; we could only set an upper bound on it (Section 5.3.4). Hence, in this Appendix 
we will be able to set only an upper limit on s to get at least some evaluation of the 
possible range of its values. We will do it in the following way. 

Equation (1) and the exact knowledge of the coefficient k in the formula for the rate 
of primaries as a function of their energy R(E p ) = kjEt can be used to calculate 
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Core locations 




400 



Figure 18: Core locations for 1702 showers giving rise to muon triggering. The opening 
angle at the antenna vertex is 2 radians. The radii of the inner and outer circular arcs 
are 100 and 400 m, respectively. Axes relative to center of array are x (geographic East) 
and y (geographic North), in meters. 

the rate R(£ v ew > The result will certainly depend on the calibrating factor s. 

Comparing this rate with the experimental upper bound, R uv {£ v ew > £°), an upper 
limit will be placed on the calibrating factor. 

In the process of evaluating R(£ vEW > £°) we will need to integrate the rate of 
primaries R(E P ) over solid angle, energy, and detection area. The question of detection 
area is of great importance since its lack of symmetry with respect to the antenna 
compelled us to make an important simplifying assumption. 

Detection area 

Core locations of showers that fire muon triggering are shown in Fig. 18. Their dis- 
tribution over the outlined area appears to be approximately uniform. To simplify the 
integration over solid angle, area, and energy, let us consider a bigger area: part of the 
ring with inner radius of 100 m, outer radius of 400 m, inside the same angle of 2 radians. 
The rate of showers passing through this area is bigger than through the one shown in 
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Fig. 18. Integration over the bigger area will lead to a greater value for R(£ uEW > 
However, that value will differ from the one calculated after integration over the outlined 
area by less than 5%, as we show below in Calculation. 

As was mentioned in Section 4.6, systematic studies have been performed only for 
simulated pulses with 5 = 5 ns. Thus, R up (£ u ew > £°) [Fig. 15(b)] was obtained under 
the assumption R ~ 200 m. Though in reality the detection area is not limited to the 
area around 200 m, we are going to accept that result as a reasonable approximation. 

Unfortunately, the area under consideration is not the full ring. If it were, it could 
have been proven that the angular distribution of the electric vector in the horizontal 
plane is approximately uniform. Then it would have been easy to derive R{E v ew > 
£°) after calculating R(£ v > In that case we wouldn't have to worry about the 
polarization of radiation from individual showers. 

In our case, however, the area is not symmetric. The general approach of calculating 
£vEW for showers with different zenith angles, azimuth angles and core locations and 
then integrating would be too difficult to implement. Therefore, it would be reasonable 
to simplify the problem by making the assumption that all showers are vertical. This 
will lead to a rough estimate of the upper limit that can be placed on the calibrating 
factor. 

Now we are almost ready to perform integration. Only one thing is missing: Since 
R up (£vEW > £„) was calculated for the EW projection of the electric field vector, we will 
need a relation between £ u and £ u ew f° r showers with different core locations. For this 
purpose let us turn to the question of polarization of radiation produced by a shower. 



Polarization 

The conventional approach assumes linear polarization in the direction perpendicular to 
both shower axis and magnetic field vector. The source of such radiation is the transverse 
current of shower particles jTU]. Alternatively, one can consider radio emission due to 
particles' acceleration in the Earth's magnetic field. For the time being we will pursue the 
latter approach but at the end of this Appendix we will give results for both alternatives. 

The magnitude and direction of the electric field vector for a radiating particle are 
governed by the general formula 



E(x,t) 



n 







7 2 (l-/3-n) 3 Z 2 



ret 



e 
c 



n x (n — f3) x /3 



1 -/3-n) 3 / 



(Al) 



ret 



where (3 is the velocity vector in the units of c, j3 = d{3/dt is the acceleration vector, 
divided by c, n is a unit vector from the radiating particle to the antenna, and I is the 
distance to the particle jlU]. The square brackets with subscript "ret" mean that the 
quantities in the brackets are evaluated at the retarded time. (A.l) does not consider 
the influence of the index of refraction of air. Here we are primarily interested in the 
direction of the electric field vector and for R > 100 m the effect of the refraction index 
on polarization is insignificant. Ref. |3Sj will give a full calculation of shower radiation, 
employing a modified formula (A.l) to take into account the refraction index. 
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Figure 19: Geometry of a vertical shower. Axes relative to antenna are x (magnetic 
West), y (magnetic South) and z (up). Vector B lies in the yOz plane. 



The first term in (A.l) decreases with distance as l/l 2 and represents a boosted 
Coulomb field. It does not produce any radiation. The magnitudes of two terms in 
(A.l) are related as l/(7 2 /) and \/3\/c. The characteristic acceleration of a 30 MeV 
electron (7 ps 60) of an air shower in the Earth's magnetic field (B ps 0.5 Gauss) is 
|a| = ecB/{^m) ~ 4.4 • 10 13 m/s 2 . Even when an electron is as close to the antenna 
as 1000 m, the first term is three orders of magnitude smaller than the second and can 
be neglected. The second term falls as l/l and is associated with a radiation field. It 
describes the electric field of a single radiating particle for most geometries relevant to 
extensive air showers. 

Consider the frame centered at the antenna, with axis Ox going to the magnetic West, 
Oy to the South and Oz directly up. For vertical showers in this frame (3 = (0, 0, — 1), 
while (3 is parallel to Ox, or, in other words, to the (1, 0, 0) vector (see Fig. 19). Let ip be 
the angle between Ox and the direction to the shower core, R the distance to the core, 
and h the altitude of the radiating particle. Then n = (~ ^g^p^ f? — J^ 1 "]^ , — ^ h ^ +R ^ j ■ 
Typical values for R are a few hundred meters, and for h several kilometers, so R/h can 
be considered small. The denominator of the second term of Eq. (A.l) is independent of 
tp. The numerator determines that, to leading (second) order in R/h, the electric field 
vector lies in the horizontal plane and is parallel to (cos 2ijj, sin 2ijj, 0). The magnitude 
of the numerator is independent of the angle ip up to terms of the order R 4 /h 4 . 
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This result shows that although particles are accelerated by the Earth's magnetic 
field in the EW direction, the polarization of the resulting radiation does not show 
preference for a particular direction. In other words, assuming uniform distribution of 
shower cores around the antenna (uniform distribution of ip), one obtains a uniform 
angular distribution of the electric field vector. 

The most important result of this section for the following calculation is the relation 
between £ v ew an d £ v '- 

£uEW = £v\ COS 2lp I 



Calculation 

Let us now calculate the rate R{£ v ew > £°) as a function of s. For vertical showers 
Eq. (1) becomes 

S v = s cos 7d exp (-^y) I* m _1 MHz _1 

where s is the calibrating factor, and j d is a dip angle (the angle between the magnetic 
field vector and the direction to magnetic North). R (is,9) only depends slightly on 9 
when 9 < 35°. The dependence on v is not very significant, either. Ref. [2] quotes two 
numbers: At v = 55 MHz R$ equals 110 m, while at v = 32 MHz R equals 140 m. The 
form of the dependence is not well known, although it is clear that Rq becomes larger as 
v decreases. Assuming linear dependence for simplicity, one can evaluate Rq = 130 m 
at v = 39 MHz, the median frequency of the main region of investigation [Section 5.3.1 
(b)]. This value of Rq will be used in this Appendix. 

Consider showers with impact parameter R, direction to the core given by angle ip, 
and primary energy 

E r >E> ( elRA)=£° 1017 eV 
v s cos 7d I cos 2^| 

where £° is expressed in jiV m _1 MHz _1 . Only these showers induce £ v ew — £v I cos 2tp\ > 
£®. So, the rate of events where the East- West projection of the pulse is greater than 
some value, R(£ u ew > £„), equals the rate Z of such showers passing through the 
detection 

Now let us determine the limits of integration. We assume that although we consider 
vertical showers only, the solid angle of observation is limited by a zenith angle 9 m = 50°. 

We also take into account that magnetic North at Dugway is located 14° east of 
true North. Therefore, the detection area is bound by lines ip = 57° + 14° = 71° and 
jj, = -57° + 14° = -43°. 

We know that the rate of primaries (in eV" 1 km~ 2 sr _1 h _1 ) is given by R(E P ) = k/Ep. 
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Hence, the integration over solid angle, area and energy gives the desired rate Z (in h 1 ) 

+00 

Z = JdflJdA J R(E P ) dE p = 

9 m 71° 400 m +00 

2tt J sm9d9 J dip J RdR J jk dE p = 

-43° 100 m E p (£°,R,4>) p 

71° 400 m 
27r(l-cosfl m ) J dip J RdRj^ = 
-43° 100 m p 
. s 2 71° 400 m 

7r£;(l-cosfl m ) / cos 2 2iP dip J Re~ m ' Ro dR (A2) 

Both integrals in Eq. (A. 2) can be readily calculated to obtain the following expression: 

R(£ uEW > £° v ) = Z^ nk(l - cose m ) (^|) ■ 2.0 • 10 3 h" 1 

where k is in eV 2 m _2 sr _1 h _1 and £® is in /iV/m/MHz. Let us note here that integration 
over the area outlined in Fig. 18 (instead of the part of the ring) has also been calculated 
numerically with the help of Mathematica, with the result 



R(£vEW > £°) « tt/c(1 - cos0 m ) ^ 



5COS7d \ 2 

10 17 £,° ) 



■1.9-10 3 h" 1 



i.e., only a 5% difference. Thus, the rapidly decreasing exponential suppresses the con- 
tribution at large distances from the antenna where the part of the ring does not overlap 
the outlined area. 

Now we can place an upper limit on s using this formula for the rate R{E v ew > £°) 
and the experimental upper limit on this rate, R up (£ uE w > £°) from Section 5.3.4: 

R{E vEW > Si) « 7rfc(l - cos6 m ) (^^) ■ 1.9 ■ 10 3 < R up {E vEW > £° v ) (A3) 



10 17 £° 



10 17 £,° 

s < 



R up (^ew > £°) (A4) 



cos7rf \ nk(l — cos6* m ) • 1.9 • 10 3 



/From Eq. ©, £°- 'R u p{£ uEW > W) ~ 0-943 /iV/m/MHz/h 1/2 . This value, substituted 
into (A.4) together with k = 3.84- 10 24 eVV^sr-V 1 = 1.38 • 10 28 eV 2 m- 2 sr~ 1 h" 1 , 
9 m = 50° and 7^ = 68°, gives the final result: s < 46. A similar calculation employing 
North-South results [Eq. (0)] gives s < 34. 

If the main source of radiation is transverse current, not particles' acceleration, then 
E is parallel to {3 x B [TJ1] and for vertical showers £ vE w = E v . This changes the integral 
of cos 2 2ip over ip in formula (A. 2) to the integral of 1 over ip, leading to s < 31. This is 
the smallest upper limit we can place on s. Unfortunately, this upper limit is not quite 
small enough to rule out the initial Haverah Park group claim of s = 20, let alone the 
values of 9.2 or 1.6. To facilitate comparison of results of different studies, let's rewrite 
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formula (A. 3) under the assumption of the transverse current mechanism. In simplified 
form it is 

> « 9.16 • 10" 4 (j^j 

The transverse current mechanism cannot account for the NS component for vertical 
showers so we will employ the model discussed above in Polarization. According to it, 
£uNS = £J\ sin 2^ | . Then, formula (A. 3) transforms into 

R(£ uNS > El) » 4.86 • 10- 4 (j^j 

The plots of R(£ v ew > £°) and R(£ uNS > £°) for different values of s are shown in 
Figs. 16 and 17. 

The noise level at Dugway is too high and the acquired sample is too limited in 
statistics and dynamic range to allow us to place upper limits strict enough to check 
the claims of the two groups. We hope that further improvement of the data processing 
technique will reduce the noise contribution. The dip angle 7^ is much smaller at the 
Auger site in Argentina (34° versus 68° at Utah). This leads to bigger electric fields for 
vertical showers, facilitating the detection of shower radiation. We also expect that the 
Argentina site will be quieter than Dugway and some clarity as regards the calibrating 
factor will be established. 



Appendix B: Relation between simulated pulse strength 
and peak detected voltage 

Within the bandwidth Sis, the antenna feels the electric field £ v • 5u. Then, the instan- 
taneous power flow is given by [£ v • 5z/) 2 /(1207r), where 1207T fl is the impedance of free 
space. The antenna effective aperture is equal to G ant (A 2 /47r), G ant ~ 2.5 being the 
antenna gain. Hence, it supplies a power of 

^-7120^ ^ 

to the input of the feedline cable. Subsequently, this power is transmitted by 60 feet 
of RG-58U cable with average attenuation at 39 MHz estimated to be 1.44 dB or a 
factor of L ~ 1.4 in power. So, the power at the input of the preamplifier is W/L. 
The impedances in the antenna system were matched, in particular, the resistance of 
the preamplifier was equal to the characteristic impedance of the cable, R = 50 Q. The 
power at the input of the preamplifier is transformed into a peak voltage V pk across the 
preamplifier resistance R, i.e. W/L = V^JR. Thus, 



. 480tt 2 v 2V pk / 120tt 2 L _ u 2V pk 120n 2 

L " r ' V G ant X 2 bv c V G ant R bv c V GR 
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where G = G ant /L ~ 1.8 is the antenna-cable system gain. For specific values used in 
data analysis (Sec. 5.3.1.) bv = 54 - 24 = 30 MHz, v = 39 MHz, 

£„ = 0.042 ¥f= = 0.031 Vpk ^V/m/MHz 
VG V / G7L8 

where V p k is in /iV. 

The above derivation implies that the electric field perceived by the antenna is due to 
frequencies in the 5u range only. However, the simulated pulse voltage with V p k = 1.3 mV 
has components in the whole frequency range (see Fig. 4, bottom), i.e., it is the result of 
an imaginary "electric field" with very diverse frequencies. To find what V p k would be 
if our antenna were sensitive to 24-54 MHz only, the forward and then inverse (for 24- 
54 MHz region) Fourier transforms of the simulated pulse voltage have been performed. 
The contribution of 24-54 MHz region into the pulse peak voltage of 1.3 mV appeared 
to be equal to 442 /iV. This leads to a corrected formula: 

£„ = Q.Q31 -j^j- Vpk =0.011-^= /iV/m/MHz 
1300 V / G/L8 ' ' 

This formula establishes the relation between the simulated pulse strength £ v and 
its peak voltage V p k at the filter-preamplifier configuration. G = 1.8 was assumed 
throughout the paper. 



Appendix C: Cost considerations for Auger project 

At present we can only present a rough sketch of criteria for detection in the 30-100 MHz 
range. Data would be digitized at a 500 MHz rate at each station and stored in a rolling 
manner, with at least 20 microseconds of data in the pipeline at any moment. Upon 
receipt of a trigger signaling the presence of a "large" shower (> 10 18 eV), these data 
would be merged into the rest of the data stream at each station. 

Per station, we estimate the following additional costs, in US dollars, for RF pulse 
detection: 



Two antennas and protection circuitry: 


200 


(a) 


Mounting hardware: 


100 


(b) 


Cables and connectors: 


200 


(c) 


Preamps: 


500 


(d) 


Digitization and memory electronics: 


2000 


(e) 


Total per station: 


3000 


(f) 



(a) Two military-surplus log-periodic antennas; crossed polarizations. 

(b) Highly dependent on other installations at site. Antennas are to be pointed vertically 
but optimum elevation not yet determined. 

(c) Antennas are mounted near central data acquisition site of each station, but suffi- 
ciently far from any sources of RF interference such as switching power supplies. Alter- 
native location of receiving stations at interstitial positions in the array would require 
microwave communication links and auxiliary power supplies and would add to cost. 
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(d) Commercial GaAsFET preamps and gas discharge tubes. 

(e) Subject to prototype development experience. Power requirements not yet known. 

(f) The number of stations to be equipped with RF detection will depend on further 
prototype experience. 

The above estimate assumes that one can power the preamps and DAQ electronics 
from the supply at each station without substantial added cost. It also assumes that 
a "large-event trigger" will be available at each station. One consideration may be the 
acquisition of antennas robust enough to withstand extreme weather (particularly wind) 
conditions. 

For detection at frequencies above or below 30-100 MHz, the criteria are not yet well 
enough developed to permit any cost estimate. 
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